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Abstract. The atomistic to continuum interface for quasicontinuum energies exhibits nonzero 
forces under uniform strain that have been called ghost forces. In this paper, we prove for a 
linearization of a one-dimensional quasicontinuum energy around a uniform strain that the effect of 
the ghost forces on the displacement nearly cancels and has a small effect on the error away from the 
interface. We give optimal order error estimates that show that the quasicontinuum displacement 
converges to the atomistic displacement at the rate O(h) in the discrete 1°° and w 1 ' 1 norms where h 
is the interatomic spacing. We also give a proof that the error in the displacement gradient decays 
away from the interface to O(h) at distance 0(h\ \ogh\) in the atomistic region and distance O(h) 
in the continuum region. E, Ming, and Yang previously gave a counterexample to convergence in 
the w 1 ' 00 norm for a harmonic interatomic potential. Our work gives an explicit and simplified 
form for the decay of the effect of the atomistic to continuum coupling error in terms of a general 
underlying interatomic potential and gives the estimates described above in the discrete £°° and 
w l,v norms. 



1. Introduction 

The quasicontinuum method (QC) reduces the computational complexity of atomistic simulations 
by replacing smoothly varying regions of the material with a continuum approximation derived from 
the atomistic model [5-9,11-15,17,18,21,23,25]. This is extremely effective in simulations involving 
defects, which have singularities in the deformation gradient. In such simulations, a few localized 
regions require the accuracy and high computational expense of atomistic scale resolution, but 
the rest of the material has a slowly varying deformation gradient which can be more efficiently 
computed using the continuum approximation without loss of the desired accuracy. Adaptive 
algorithms have been developed for QC to determine which regions require the accuracy of atomistic 
modeling and how to coarsen the finite element mesh in the continuum region [1-3,16-18,20]. The 
atoms retained in the atomistic region and the atoms at nodes of the piecewise linear finite element 
mesh in the continuum region are collectively denoted as representative atoms. 

Recent years have seen the development of many QC approximations that differ in how they 
compute interactions among the representative atoms. In the following, we concern ourselves with 
the original energy-based quasicontinuum (QCE) approximation [15,25], but the phenomena that 
we analyze occur in all other quasicontinuum approximations, as well as in other multiphysics 
coupling methods [8]. In QCE, a total energy is created by summing energy contributions from 
each representative atom in the atomistic region and from each element in the continuum region, 
where the volume of the elements in the atomistic to continuum interface is modified to exactly 
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conserve mass. This construction was chosen so that for any uniform strain the QCE energy, the 
continuum energy, and the atomistic energy are identical. (As discussed later, this conservation 
property for the QCE approximation is not sufficient to prevent the existence of nonzero forces at 
the atomistic to continuum interface for uniform strain.) The representative atoms then interact 
via forces defined by the total energy. This makes for a simple and versatile method that can 
treat complicated geometries and can be used with adaptive algorithms that modify the mesh and 
atomistic regions during a quasi-static process. Other atomistic to continuum approaches have 
been proposed, for example, that utilize overlapping or blended domains [4,19]. 

One drawback of the energy-based quasicontinuum approximation that has received much atten- 
tion is the fact that at the atomistic to continuum interface the balance of force equations do not 
give a consistent scheme [22]. As explained in Section [21 the equilibrium equations in the interior 
of both the atomistic region and the continuum region give consistent finite difference schemes for 
the continuum limit, whereas the QC equilibrium equations near the interface are not consistent 
with the continuum limit. This is most easily seen by considering a uniform strain, which will be 
assigned identically zero elastic forces by any consistent scheme. (Ensuring that a given scheme 
assigns zero forces for uniform strain has been known as the "patch test" in the theory of finite 
elements [24].) The nonzero residual forces present in QCE for uniform strain have been called 
"ghost forces" [7,22]. 

In this paper, we give optimal order error estimates for the effect of the inconsistency on the 
displacement and displacement gradient for a linearization of a one-dimensional atomistic energy 
and its quasicontinuum approximation. We consider the linearization of general interatomic po- 
tentials which are concave near second-neighbor interatomic distances. This property guarantees 
that the interfacial error due to the Cauchy-Born approximation with a second-neighbor cut-off is 
positive [7, p. 117] and that the quasicontinuum error is not oscillatory in the atomistic region (see 
Section [3]). Similar optimal order error estimates have been given by E, Ming, and Yang [10] for a 
harmonic interatomic potential. 

We begin by linearizing a one-dimensional atomistic energy, its local quasicontinuum approxi- 
mation (which we will call the continuum energy), and its quasicontinuum approximation about a 
uniform strain for a second-neighbor atomistic energy. We will show in Section [2] that the three 
systems of equilibrium equations are then 

f, (atomistic) 
f, (continuum) 
f, (quasicontinuum) 

where f is an external loading, L and u are the linearized operator and corresponding displacement 
for each scheme, g is non-zero only in the atomistic to continuum interface, and h is the interatomic 
spacing. The term g in the quasicontinuum equilibrium equations is due to the unbalanced second- 
neighbor interactions in the interface f)2. 13j) and for uniform stretches is precisely the ghost force 
described in [7,15,22]. 

Formally, the error decomposes as 

u a - u 9C = {(L a,h )~ l - (L^y^f - (L gc ' ,t )~ 1 g. 

(The operators are all translation invariant, so they only have solutions up to the choice of an 
additive constant.) In this paper, we focus on the second term, (L (?c ) _1 g, which is the error due to 
the inconsistency at the interface. To do so, we consider the case of no external field, f = 0, which 
will make u a = 0. For most applications of the quasicontinuum method, the only external field is 



L a > h u a 



QC 



g 
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due to loads that are applied on the boundary of the material, far from the atomistic to continuum 
interface. 

We showed in [7] that the ghost forces are oscillatory and sum to zero. In this paper, we prove that 
the error in the displacement gradient is 0(1) at the interface and decays away from the interface 
to 0(h) at distance 0(h\ \ogh\) in the atomistic region and distance 0(h) in the continuum region. 
As noted above, similar results have been given in [10] for a harmonic interatomic potential with 
f 7^ and Dirichlet boundary conditions. Here, we present a simplified approach starting from a 
linearization of a quasicontinuum approximation with a concave second-neighbor interaction. We 
explicitly give the form of the solution and analyze the solution in discrete l°° and w 1,p norms. We 
show that the quasicontinuum displacement converges to the atomistic displacement at the rate 
0(h) in the discrete l°° and w 1,1 norms where h is the interatomic spacing. 

In Section [2j we describe the energy-based quasicontinuum approximation (QCE) and set up 
the analysis. In Section [3j we prove Theorem 13.11 for the quasicontinuum energy that gives an 
optimal order, 0(h) error estimate in the l°° norm and a 0(h l l p ) error estimate in the w l,p norm 
for 1 < p < oo. Note that for simplicity the models and analysis are presented for the case where no 
degrees of freedom have been removed in the continuum region, but we explain in Remark 13. II that 
identical results hold when the continuum region is coarsened. We present numerical computations 
in Figure Q] that clearly show that the error is localized in the atomistic to continuum interface. 



2. One-Dimensional, Linear Quasicontinuum Approximation 

We consider an infinite one-dimensional chain of atoms with periodicity 2F in the deformed 
configuration. Let yj denote the atomic positions for — oo < j < oo, where there are 2N atoms in 
each period. Let h = 1/N and let 

u j '■= Vi ~ Fh 3 

denote the displacement from the average interatomic spacing, Fh. In the following, we analyze 
the behavior of the quasicontinuum method as the atomistic chain approaches the continuum limit 
with F fixed and N — > oo. 

The atomistic energy for a period of the chain is 



i=-N+i L V ' / V " / 



fjVj 



(2.1) 



where <p(r) is a two-body interatomic potential (for example, the Lennard-Jones potential <p(r) = 
1/r 12 — 2/r 6 ) and f = (f-N+i, ■ • • , In) are external forces applied as dead loads on the atoms. The 
periodic conditions 

Vj+2N = Vj + 2F or u j+2 N = Uj 

allow £ tot ' h to be written in terms of y := (y_/\r+l> • • • iVn)- We assume that fj = ^ 

otherwise there are no energy minimizing solutions since the elastic energy is translation invariant. 
In the following, we discuss the existence and uniqueness of solutions to each of the models we 
encounter. We note that the energy per bond in (|2.1j) has been scaled like hcp(r/h). This scaling 
implies that if we let yj = y(j/N) and fj = f(j/N) for j = -N + 1, . . . , N where y G C^Q-l, 1]) 
and / £ C([— 1, 1]), then as N — > oo and F is held fixed, the energy of a period (|2.ip converges to 

I 

Hy'( x )) + <K2y'(z)) - f(x)y(x) dx. 

l 
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We expand first neighbor terms around F, giving 
, ( Vj+i ~ Vj ' 



h 



F+ u j+ t 



h 



<j>(F) + <j>'(F) 



h 



u j+ i 



+ o 



Uj+l 



and the second neighbor terms around 2F, giving 



^V. - fit 



h 



2F + 



u j+2 



<P(2F) + <P'(2F) 



Un 



h 



+ 



1 Jjl 



(2F) 



h 



+ o 



U j+2 



2.1. Atomistic Model. The linearized atomistic energy is then given by 

2 



£ a ' h (u) 



h 



N 

£ 

j=-N+l 



h 



J -4- 1 rh" 

r 



U j+1 



u 



U j+2 - Uj i „ 
T r 2 ( r2F 



h 



h 



fjUj 



(2.2) 



where </>' F := <j>' (F) , ^ := <f/'(F), <f>' 2F := </f(2F),<% F := <f/'(2F), and u := (u- N+1 , . . . ,u N ). Note 



that here and in the following, we neglect the additive constant <ft(F) + (j>(2F) — hJ2f=-N+i fjFhj 
in the linearized energy. We assume that <fi E C 2 ([ro, oo)) for some tq such that < ro < F, and 

<jJ' F > and <j%p < 0. (2.3) 

This holds true for the Lennard- Jones potential for Fh below the load limit, unless the chain is 
extremely compressed (less than 60% of the equilibrium length). The property (f>2 F < ensures 
that the quasicontinuum error is not oscillatory in the atomistic region (see Section [3|) . 
We furthermore assume that 

<jt' F + 5$p > 0, (2.4) 
which will be sufficient to give solutions to the QC equilibrium equations under the assumption of 
no resultant external forces (see Lemma |2.1|) . In contrast, the weaker assumption <p F + 4:<fi'n F > 
is sufficient for the fully atomistic or fully continuum approximation. The equilibrium equations, 

l d£ c ' 



h du 



— (u) = 0, for the atomistic model (12. 21} are 



(L^U); 



-^2 F Uj +2 ~ <j>pUj + l + 2(4>"p + (p2p)Uj - (tfpUj-X - (f) 2 'pUj-2 



(2.5) 



h 2 

Uj+2N = Uj, 

for — oo < j < oo. Note that scaling by j- makes this a consistent approximation of the boundary 
value problem 

— (<j>'p + 4(j)2p)u" (x) = f for — oo < x < oo, 
u(x + 2) = u(x) for — oo < x < oo. 

The linearized atomistic energy ()2.2[) has a unique minimum (up to a constant) if (j)" F + 4</>^p > 0, 
provided that Ylf=N-l h = Standard ODE results show that (|2.6|) has a unique solution (up to 
a constant) provided that J_ x f(x) dx = 0. 



(2.6) 
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Remark 2.1. For the atomistic energy (|2.2p . the linear terms sum to zero by the periodicity of the 
displacement, since 



N 

» £ 

j=-N+l 



» F - 



+ 6' Uj+2 
+ 92F ' 



h ' h 

F [UN+1 ~ U-N+l] + <t>2F [uN+2 + UN+1 ~ ~ U_jV+l] 



However, we keep these terms in the model since they do not sum to zero when the atomistic model 
is coupled to the continuum approximation in the quasicontinuum energy. The resulting terms give 
a more accurate representation of what happens in the non-linear quasicontinuum model. 

2.2. Continuum Approximation. The continuum approximation splits the chain into linear 
finite elements with nodes given by the representative atoms, which we recall are a subset of the 
atoms in the chain. The energy of the chain is the sum of element energies which depend only 
on the element's deformation gradient, the linear deformation that interpolates its nodal positions. 
The energy of an element is then computed by applying the element's deformation gradient to the 
reference lattice, computing the energy per atom using the atomistic model, and multiplying by the 
number of atoms in the element (where the boundary atoms are shared equally between neighboring 
elements). If the continuum approximation is not coarsened (every atom is a representative atom), 
then the continuum energy is given by 



£ c ' h (u) 



N 

■ £ 

i=-iv+i 



+ 202 



2F I 



U j + 1 



Un 



+ 



+ m 



2F) 



3+1 



h u 3 



(2.7) 



See [7] for a derivation of the continuum energy and a discussion of the error terms at the element 
boundaries. For j € {— N + 1, . . . , N}, the equilibrium equations for the continuum approximation 
are 



U 



f>F + ^2F) 



i + 2uj 



ij-i 



h 2 



fi, 



(2- 



which is also a consistent approximation for the boundary value problem (|2,6p . It is easy to see that 
the continuum energy (|2.7p has a unique minimum (up to a constant) if <p" F + 40op > 0, provided 
that Ylf=N-i fj = 0- The quasicontinuum method inherently supports coarsening, but we neglect 
it here since in one dimension this only changes the scaling of equilibrium equations. 



2.3. Splitting the Energy. We can split the atomistic energy and the continuum energy into 
per-atom contributions so that 



£ 



a,h i 



U 



N 

Yl [ £ j' h ( u ) - fi u i 

=-N+l 



and £ c 



h. 



u 



N 

E 

=-N+l 



fi u 3 



There are many possible ways to define the per-atom contributions, and we do this in such a way 
that these contributions are linearizations of the ones in the fully nonlinear case presented in [7, 25] . 
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In this case, we split the energy of each bond to obtain 



6 



c a,hr \ 1 



u i+i ~ u i ,i in / u j+i 



u. 



+ Wf 



, jj u j+i - U J L i jj' ( u i+ 2 ~ u i 

' 92F T~ 



1 

+ 2 



+ 2V2F 



f Uj Uj-1 i ,n ( Uj Uj-l 
T) F h 



, Uj - Uj-2 i ,„ / - «j-2 
+ 7 r 2 < +'2F 



(2.9) 



and 



£?*(u) := 



(0' F + 2<£' 2F ) 



+ 



&' F + 2<j)' 2F ) 



Uj - u,_i 



(2.10) 



2.4. Energy-Based Quasicontinuum Approximation. The energy-based quasicontinuum ap- 
proximation partitions the representative atoms into atomistic and continuum representative atoms 
and assigns to each atom the split energy corresponding to its type (j2.9M2.10"]) . We define the nodes 
— N + 1, ... j —K — 1 and K + 1, . . . , N to be continuum and —K, . . . , K to be atomistic, where 
we assume that 2 < K < N — 2 to ensure well-defined atomistic and continuum regions. The 
quasicontinuum energy is then 



-K-i K N 

£^ h (u):= £ (u)+ £ £f(u)+ Y, £ f 2J fj^r 

j=-N+l j=-K j=K+l j=~N+l 



(2.11) 



Since the energy is quadratic, the equilibrium equations, ^ ^f^~ ( u gc) = 0; take the form 



L^ h u qc - g = f . 



(2.12) 
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For < j < N, the QCE operator is given by 
{L qc ' h u)j 



+ < 



Jl ~ u j+l + 2u j - u j- 


i 




up 


7 9 






\^2F 


—Uj+2 + *Uj — 


Uj 


-2 




Ah 2 






Acp 2F 


—Uj + 2 + 2Uj — 


Uj 


-2 


Ah 2 






H2F 


—Uj + 2 + 2Uj - 


Uj 


-2 


Ah 2 






±<t>2F 


—uj + i + 2uj — 


Uj 


-1 


h 2 








—uj + i + 2uj — 


Uj 


-1 




h 2 






±<f>2F 

V 


—Uj+i + 2uj — 


Uj 


-1 


h 2 







+ 



<P 2F Uj +2 ~ Uj 








h 2h 








2(f>2 F Uj+l - Uj 


+ 




Uj +2 ~ Uj 


h h 


h 


2h 


2(f) 2 p Uj — Uj-i 


+ 


1>2F 


Uj - Uj-2 


h h 


h 


2h 


<f>2p Uj - Uj-2 








h 2h ' 









0<j<K-2, 
K 
K, 
K 

j = K + 2, 

K + 3 < j < N. 



1, 



1, 



Similarly, g is given by 



9j 



0, 

-\(j) 2F /h, 
\<t>' 2F /h, 

\4>' 2F /K 
-\cj)' 2F /h, 

0, 



< j < K - 2, 
j = K-l, 
j = K, 
j = K + l, 
j = K + 2, 
K + 3 < j < N. 



(2.13) 



For space reasons, we only list the entries for < j < N. The equations for all other j £ Z 
follow from symmetry and periodicity. Due to the symmetry in the definition of the atomistic and 
continuum regions, we have that L qc - h = L q ^_- and gj = —g~j for —N + 1 < i, j < 0. To see this, 
we define the involution operator (Su)j = —U—j and observe that £ qc,h (Su) = £ qc,h (u). It then 
follows from the chain rule that 



S T L qc,h Su _ S T g _ S T { = L qc,h u _ ( 

bmce we can conclude that 

SL qc > h S = L qc ' h and Sg = g. 



f for all periodic u and f . 



(2.14) 



Note that the expression for g does not depend on (p' F since the first-neighbor terms identically 
sum to zero in the energy (12. lip . We can now observe that the QCE approximation (|2.12p is not 
consistent with the continuum limit of the atomistic model (|2.6p . 

The linear operator L qc has all uniform translations, u = cl = (c, c, ...,c), in its nullspace. 
To see that this is the full nullspace, we consider the factored operator L qc = D T E qc D, where 
(Du)j = ^±pi and 



(E qc r)j 



<P2F r 3-l + 
4>2F r 3 
2<f>2F' J 
2<f>2F r j 



iff 



i + 

Tj-l + 



F + 
n 
F 



{ 2F>3 + 



' ) 2F r j+^ 



Ul _i 3 ah 



2^2F) r i 



+ \$2F r 3+^ 

rj + l^ F r j+1 , 



Ul 



+ ■ 



-1 + 

Iff • 



2F) r V 



< j < K - 2, 
j = K-l, 
j = K, 
j = K + l, 
K + 2 < j < N. 
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We see that E qc is diagonally dominant provided <j)' F + 5<f>2 F > 0, hence assumption (|2.4p implies E qc 
is invertible. So we have that the nullspace of L qc is precisely the nullspace of D. Thus, L qc u = g 
has a solution whenever Ylj=-N+l fj = 0> since Ylj=-N+l 9j = ®- This solution is unique up to a 
constant. 

We now gather together the existence and uniqueness results stated for the models. 

Lemma 2.1. If Ylj=-N+i fj = an d 4>"f + > 0, ^ en ^ e linearized atomistic energy ()2.2[) 

and continuum approximation ()2.7[) 6o£/i /iaue a global minimum that is unique up to an additive 
constant. 

Under the slightly stronger assumption (j)" F + hcjJ^p > 0, t/ie quasicontinuum energy ()2.1ip /ias a 
unique minimizer up to a constant. 

Here, and in the following, we take f = 0, in order to focus on the effect of the ghost force g. 
Under this assumption, we can conclude that the unique mean zero solution to the QCE equilibrium 
equations (12. 12H is odd. This follows from = S and (I2.14p which together imply that Su is a 
solution if and only if u is. Because S preserves the mean zero property, we conclude that u qc is 
odd. The unique odd solution to the atomistic equations, L a ' h u a = 0, is u a = 0. Thus, the QCE 
equilibrium equations, 

L* C >V - g = 0, ( 2 - 15 ) 
are also the error equations, and the quasicontinuum solution is the error in approximating u a . 

2.5. Discrete Sobolev Norms. The effect of the interface terms on the total error is norm- 
dependent, so we now employ discrete analogs of Sobolev norms [18]. We define the discrete weak 
derivative by 

, Uj + l - Uj 

For 1 < p < oo the discrete Sobolev norms are given by 



h 




and for p = oo by 



u p = max \u 
h -N+l<j<N 

MLi-°° = \\v\\ts° + 



3h 



"h 

The above discrete Sobolev norms are equivalent to the standard Sobolev norms restricted to the 
continuous, piecewise linear interpolants u(x) satisfying u(j/N) = Uj for j = —N + 1, . . . , JV. 

3. Convergence of the Quasicontinuum Solution 

We now analyze the quasicontinuum error, u qc . We note that is it theoretically possible to 
solve (|2.15p explicitly for u qc ; however, the form of the solution is complicated by the second- 
neighbor coupling in the atomistic region, so we instead obtain estimates for the decay of the 
error, u gc , by analyzing a 0(/i 2 )-accurate approximation of the error. Figure [1] shows the results of 
solving (|2.15p numerically for odd solutions, Uj = —U—j, with three choices of lattice spacing and 
two sets of parameters. Note that for both sets of parameters, the magnitude decays linearly with 
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h, whereas the displacement gradient is 0(1) in the atomistic to continuum region. The following 
argument proves the qualitative error behavior analytically. 

3.1. Form of the Solution. In the interior of the continuum region the solution is linear, but in 
the atomistic region u qc is the sum of a linear solution and exponential solutions. The homogeneous 
atomistic difference scheme 

— 4>2 F Uj + 2 - <j)'pUj+i + (2<j)"p + 24>2 F )Uj - <p" F Uj-i - <p2 F Uj-2 = (3.1) 

has characteristic equation 

-^ F A 2 - <^A + (20'; + 2$ F ) - ^A" 1 - ^ F A~ 2 = 0, 



with roots 



where 



1) 1) A, 



A 



-2</>» F 

Based on the assumptions on in f|2.3j) and (|2.4p . we have that A > 1. We note that if §'L F were 
positive contrary to assumption (12. 3D . then A would be negative which would give a damped oscilla- 
tory error in the atomistic region. General solutions of the homogeneous atomistic equations (13.1j) 
have the form Uj = C\ + C^hj + C^X 3 + C^A - - 7 , but seeking an odd solution reduces this to the 
form Uj = C 2 hj + C7 3 (A- 7 ' - A~ J ). 

The odd solution of the quasicontinuum error equations (|2. 15[) is thus of the form 

'mxhj + pi^P-), 0<j<K, 
( u qc)j = \ rri2hj — 771,2 + UK+1, j = K + 1, 

,7772^- 7772, K + 2 < j < N, 

where expressing the unknown uk+i using a perturbation of the linear solution, uk+i, simplifies 
the solution of the equilibrium equations. The four coefficients mi, ni2, uk+i, an d P can be found 
by satisfying the four equations in the interface, j = K — 1, . . . , K + 2. Summing the equilibrium 
equations across the interface gives 



K+2 K+2 

qcjj 



= E si = E ( Lqc,hu 

j=K-l j=K-X 



,„ U K -1 ~ U K -2 , . ,// UK + UK-1 ~ U K -2 ~ U K -3 

bp — p — + 4< ^ ; — 

UK+3 — Uk+2 



b F + 4(p2 F 



h 2 



(A" j_AA!< \( mi m2 

The cancellation of the exponential terms in the final equality holds because 

cj>'i F {X K - \- K ) + + 4>2f)^ K ~ X - ^ K+1 - A^ 2 + \- K+2 ) + 2 V(-A X - 3 + A"^ +3 ) = 0, 

which can be seen by summing (|3.ip with the homogeneous solution Uj = —A- 7 for j = —K + 
2,...,K — 2. Thus mi = 1 that is, the slope of the linear part does not change across the 
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Decay of Quasicontinuum Ghost Force Effect 

h = 1/20 (j)" F = 1 (j>% F = -0.05 h = 1/20 <j>" F = 1 <f>% F = -0.15 




h = 1/40 <f>" F = 1 <j>% F = -0.05 



-0.5 0.5 

h = 1/40 4>" F = 1 4>'^ F = -0.15 



1 




0.01 
0.005 


-0.005 
-0.01 



1 -0.5 0.5 1 

h = 1/80 4> F = 1 4>'^ F = -0.05 
















-0.02 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 

FIGURE 1. Error for the energy-based quasicontinuum approximation, u qc . We 
observe that the magnitude of the error is 0(h). However, the oscillation near the 
interface means that the error in the displacement gradient is 0(1) in the interfacial 
region. The average deformation gradient, F, for the right column is close to failing 
the stability condition <j)" F + > 0. In all plots K = N/2 and (f>' 2F = 1. 
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interface. Hence, the odd solution is given by 

mhj ■ !' x 

i U qc)j 



mhj 
^mhj 



m + uk+i 



0<j<K, 
j = K + l, 
K + 2 < j < N, 



(3.2) 



where the coefficients m,UK+i, and f3 can now be found by satisfying any three of the equations 
in the interface, j = K — 1, . . . , K + 2. 

3.2. Magnitude of the Solution. We focus on the equations at j = K — 1, K + 1, and K + 2 
and split the interface equations as (Ak + hB)x = hh, where 











_1 




A K = 


-4>"f~ 


h 




2<P" F + f 


<t>2F ~ 






h 


%F 


-4>" F - ^ 


4>2F 









" 




m 


B = 










, x = 


UK+l 




<t>2F 













■\4>2FlK 



b= 



2F 



and 7j = XJ 7k ■ We note that Ar-, -B, and b do not depend on /i directly, though Ak may have 
indirect dependence if K scales with h as in Figure [TJ Therefore, we can neglect B and conclude 
that x is 0(h) provided that Aj} exists and is bounded uniformly in K. 

Lemma 3.1. For all K satisfying 2 < K < N — 2, the matrix Ax is nonsingular and \\A K \\ < C 
where C > is independent of K. 



Proof. Applying row reductions gives the upper triangular form 

-\(t>2F (t>2FlK+l ~ \4> 2 'f1K-\ 

o -4>" F - f^V 4>2 F iK+i - \4>2 F ik - I^Vm-i 

T) K 



A 



2 < r2F 



where 

m = m? + f^ F + f(^ F ) 2 ) (2 lK+ i 



IK - iK-l) + 2^2F Wf + W2f) (iK ~ 7/f-l 



If the diagonal entries of A are non-zero, then Ak is nonsingular. The coercivity assumption 
(j) F + 5<f% F > d23D implies that - 9/2$^ < since < 0, so the first and second 

diagonal entries are non-zero. Since the second term of t/k is negative, we can use the fact that 
IK ~ IK-l < 2-/K+1 ~ IK ~ 7K- 1 to see that 



'2FJ ) ( 2 7^+l 

y F + 



i 

s/2 



IK ~ 7K-1) 

{27K+1 



IK ~ IK-l) 



> 0. 



Therefore, A K l exists for all K. Taking limits, we find 



lim t]k ~> 

K— >oo 



+ 4 



l 

V2 



(2X-1- A" 1 ) > 0, 



where we note that the elementary matrices corresponding to the row reduction operations did not 
depend on K so that lim^^oo Ak is nonsingular. The inverse of a matrix is continuous as a function 
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of the entries whenever the matrix is nonsingular. Thus, the fact that lim^- 
implies that lim^-_ 

we conclude that \\A^ || is uniformly bounded 



A^- 1 !) is finite. Since | I^^H is finite for all K and lim^-. 



Ak is nonsingular 



\AjsW is finite, 



□ 



Thus, we have that iti,uk+i, and (5 are all 0(h). We can express the derivative, u' , as 



m + 
m — 
m 
m, 



( A^+^A-J- 1 



h 



+ 



h 
h 



h 



A^ 

3 ( \ K -\- K 
~ —IF 



< j < K - 1, 

j = K, 

j = K + l, 

K + 2<j<N-l, 



where u'_j_ 1 = u'j for j = 0, . . . , N — 1. 



Theorem 3.1. Let u qc be the solution to the QC error equation (|2.15p . Then for 1 < p < oo, 
2 < K < N — 2, and h sufficiently small, the error can be bounded by 

< Ch, 



u, 



qc 



||u„ c || i, P < Ch 1 ^, 

II qc\\ w <f — J 

h 

where C > is independent of h, K, and p. 

Proof. The result for the £°° norm follows from the fact that all terms in 



are 0(h). To show 



the bound on w l,p , we first apply the triangle inequality to separate the m, Uk ^ 1 , ^, and j- terms 
which we bound using the fact that UK+i,m, and (3 are 0(h). We have 




□ 

Finally, we show that the pointwise error in the derivative, u' , decays exponentially in j to 
0(h) away from the interface in the atomistic region and decays immediately to 0(h) away from 
the interface in the continuum region. 

Lemma 3.2. There is a C > such that \(u' qc )j\ < Ch for all < j < K + and K + 2 < j < N. 
Thus, the interface has size 0(h\ log h\). 



Proof. For h sufficiently small, we have that max(m, f3) < Ch. Since u'j = m for K + 2 < j < N, in 
this region u'j < Ch. For the terms < j < K — 1 it is sufficient to show that the exponential term 
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is less than or equal to Ch. For < j < K + jS-j, we have that 



In A 



X K X K 



< X ^- K 

< Ch. □ 



Remark 3.1. In order reduce the degrees of freedom, the continuum region is coarsened in compu- 
tations using the quasicontinuum method. For simplicity, coarsening was omitted from the model 
presented in this paper, but, in fact, the results are unchanged if it is used. Conventionally, coars- 
ening only occurs away from the atomistic to continuum interface, so that no degrees of freedom 
are removed if they interact directly with the atomistic region. Since the solution uj is linear for 
K + 2 < j < N, any level of coarsening produces an identical solution. 
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